Thermal Parameters in Heavy Ion Collisions at SPS and RHIC: 

Centrality Dependence 
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The centrality dependence of thermal parameters describing hadron multiplicities and 
intermediate-mass dilepton spectra in heavy-ion collisions at SPS and RHIC is analyzed. From 
experimental hadron multiplicities we deduce evidence for strangeness saturation at high energy 
and maximum centrality. The observed dilepton spectra can be parameterized by a centrality inde- 
pendent temperature. 
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I. INTRODUCTION 

It has been shown that various observables of relativistic heavy-ion collisions can be well described by statistical- 
thermal or hydrodynamical models. In such a way, a selected subset of a large number of observables can be reproduced 
by a small number of characteristic parameters, such as temperature, density or flow velocity. It is the subject of 
the present note to pursue this idea and to analyze the centrality dependence of the thermal parameters describing 
hadron multiplicities and dilepton spectra in the intermediate-mass region. This will provide further information 
about the effects of the size of the excited strongly interacting system and help in the systematic understanding of 
the experimental data. 

The aim of the present work is twofold: (i) the previous analysis ||] of hadron multiplicities at SPS and RHIC is 
repeated and improved, e.g., by including feed-down from weak decays, and (ii) the consideration of the intcrmcdiatc- 
mass dilepton spectra at SPS is outlined in more detail, in particular by analyzing both the invariant mass and 
transverse momentum spectra. 



II. HADRON MULTIPLICITIES 



Tr 



Hadron multiplicities can be reproduced 
exp{-(iJ~/i"Qf)/T} 



by the grand-canonical partition function Zi{V,T,fi'^ 



and 



where H is the statistical operator of the system, T denotes the temperature, 

/i" and represent the chemical potentials and corresponding conserved charges respectively. In the analysis of 47r 
data, the net-zero strangeness and the baryon-to-electric charge ratio of the colliding nuclei constrain the components 
— (M-BiMSiMq); where the subscripts B,S and Q refer to baryon, strangeness, electric charge. These constraints 
have to be relaxed when considering data in a limited rapidity window, increasing the number of free parameters. 
The particle numbers are given by 
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where we include phenomenologically a strangeness saturation factor 7^ (with Si the total number of strange quarks 
in hadron species i) to account for incomplete equilibration in this sector, Ei — ^^p'^ + mf, and BW is the Breit- 
Wigner distribution (to be replaced by a 5- function for stable hadrons). The final particle numbers are Ni — Nf^"^ + 
Br^^^A'J"'" due to decays of unstable particles with branching ratios Br-'^*. Such a description can certainly be 
justified for multiplicities measured over the whole phase-space, since many dynamical effects cancel out in ratios of 
hadron yields. 
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We have analyzed two data sets: (i) NA49 An multiplicities of (tt) = ^(7r+ + t^~),K^, p, (p, and -/Vpart (taken as 
the sum over all baryons) in 6 centrality bins in the reaction Pb(158 AGeV) + Pb (it should be emphasized that 
protons are not included in our analysis |^ since there may be a spectator component in non-central collisions), and 



(ii) PHENIX mid-rapidity densities of tt , K^, and in the reaction Au -I- Au at -^s — 130 AGeV in 5 centrality 
bins [0 . In 1^ it is estimated that the probability for reconstructing protons from A decays as prompt protons is 32% 
a± Pt — 1 GeV/c. This is the probability for feed-down used in our present fits to the mid-rapidity PHENIX data. 
The results of our fits are displayed in Fig. 1. Specific features of the mesons at SPS will be discussed elsewhere [||. 

A comparison of the individual thermal parameters of both data sets is displayed in Fig. 2, which corrects Fig. 2 
in jl] (the correction of a coding error in j^] lowers the temperature and brings it in agreement with |^). Most 
remarkable is the drop of the baryo-chemical potential b and the rise of the strangeness saturation factor 7^ when 
going from ^/s = 17 AGeV to 130 AGeV. It would appear that T remains essentially constant for -y/s — 17 AGeV, 
while at 130 AGeV it rises with iVpart- The parameter /is is fairly independent of the centrality, while the strangeness 
saturation factor increases with centrality. Despite the rather limited set of analyzed hadron species, the extracted 
thermal parameters describe other hadron yields, which are at our disposal in central collisions, fairly well. 
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FIG. 1. Comparison of NA49 data (left panel, symbols, 
(lines) . 



and PHENIX data (right panel, symbols, Q) with our model 
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FIG. 2. The temperature, baryo-chemical potential and strangeness saturation factor as a function of A'^part. Solid squares 
(open circles) are for the NA49 (PHENIX) data §|| ( [§). 



III. INTERMEDIATE-MASS DILEPTONS 

As pointed out in the dilepton spectra can be analyzed in the same spirit as the hadron multiplicities and 
hadron momentum spectra, i.e., one discards any detail of the dynamics and asks only for a simple parameterization. 
As a result, one gets for the thermal dilepton spectrum |^ 

dN ho? ( M_L C0sh(r - Fcms) ' 
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where Q is the lepton pair's four-momentum, M±_ its transverse mass and Y its rapidity; Fcms denotes the fireball 
rapidity, and A'dii is a normalization factor characterizing the space-time volume of the fireball; flow effects are 
negligible for invariant mass spectra. Eq. (||) is based on the quark-hadron duality . 

In we have shown that the space-time averaged temperature parameter Tdii ~ 170 McV (i.e., a value coinciding 
with the chemical freeze-out temperature) provides a common description of the low-mass CERES data |Q and the 
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FIG. 3. Fits of the thermal contribution according to Eq. (2) with free temperature parameter and normalization factor to 
the NA50 data The data are for 7 centrality bins {Et = 19, 36, 52, 67, 80, 93, 110 GeV for bin 1 - 7). The DrelbYan yield 
and the open charm contributions are calculated as described in detail in [H,ll[. 
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intermediate-mass NA50 data at ^/s = 17 AGeV; also the WA98 photon data ||T^ are consistent with this value. 
It was, therefore, a surprise to us that our analysis of the efficiency corrected and centrality binned NA50 data jl^] 
separately with A'dii and Tdii as free parameters give as optimum fit a temperature scale in the order of 250 MeV 
, see Figs. 3 and 4.[] At this temperature the spectral shapes of thermal dileptons (Therm) and dileptons from 
semileptonic correlated decays of open charm mesons (DD) are nearly identical within the NA50 acceptance. 




'^Note that in Figs. 3 and 4 the used normalization of both the Drell-Yan yield and the open charm contribution differs by 
a factor 1/1.5 from the subsequent figures. This covers the uncertainty in fitting the large-M Drell-Yan tail and the relative 
normalization of the charm contribution. 
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To illustrate the sensitivity on the temperature parameter we exhibit in Figs. 5 and 6 a comparison with the data 
when using a fixed value of Tdii = 170 MeV. In contrast to the fits displayed in Figs. 3 and 4 the normalization 
factor now systematically increases with centrality, as expected. This could be interpreted as a hint that the complete 
disentangling of temperature parameter and normalization factor is not fully adequate.^ Rather, due to the space-time 
evolution both quantities are correlated. Because of uncertainties in the normalization and the shape of the Drell-Yan 
(DY) and charm contributions the previously found value of Tdii ~ 170 MeV [^,0 can not be excluded. 
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FIG. 5. As in Fig. H but for fixed T = 170 MeV and free parameter Vdii. 



We thank R. Rapp for pointing out this possibility. 
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As in Fig. H but for fixed Tdii = 170 MeV and free parameter A'^dii- 



In addition to the invariant mass spectra also transverse momentum spectra integrated over the intermediate- 
mass region AI = 1.6 • • • 2.5 GeV are at our disposal. At first glance, they do not better constrain the temperature 
parameter T^^ since the large-Qj^ region is dominated by the Drell-Yan yield, see Figs. 7 and 8. The 1ow-Q_l region, 
however, favors a smaller value of Tdii, while the large-Q^ region could be fine tuned by a slight change of the intrinsic 
transverse parton momentum (we use here \/Jk^ = 0.8 GeV [||). 

In contrast to the invariant mass spectra, where the thermal yield at Tdii ~ 250 MeV has the same shape as the 
charm yield, the transverse mass spectra of charm and thermal contributions look similar at Tdii ~ 170 MeV. 
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IV. SUMMARY 

The analysis of hadron multiplicities indicates a centrality independence of the chemical freeze-out temperature 
and baryo-chemical potential for ^/s = 17 AGeV, while at ^/s = 130 AGeV the baryo-chemical potential is constant 
but the temperature increases with centrality. The strangeness saturation increases with centrality for both =17 
AGeV and 130 AGeV and appears to approach unity with further increasing energy and maximum centrality. 
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FIG. 8. As Fig. 7 but for Tdii = 250 MeV. x^/d-o.f. is marginally larger than for Tdii = 170 MeV in Fig. 7. 

The analysis of the intermediate-mass dilepton spectra taken by the NA50 coUaboration is not yet conclusive. Our 
parameterization of the thermal contribution is based on quark-hadron duality where the shape and the normalization 
are not interrelated. We find an indication that the space-time averaged temperature determining the shape of dilepton 
spectra stays fairly independent of centrality, but is significantly larger than previously anticipated in 1^,0 thus being 
in agreement with values deduced in from photon data. It should be noted, however, that the analysis is sensitive 
to details of the Drell-Yan and open charm dilepton yields, both their shape and normalization. In view of these 
inherent uncertainties, the centrality independent value of Tdn « 170 MeV appears still compatible with the NA50 
data. 
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